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Abstract 

An imitation of 2d field theory is formulated by means of a model 
on the hierarchic tree (with branching number close to one) with the 
same potential and the free correlators identical to those of 2d ones. 
Such a model possesses some features of original models for certain 
scale invariant theories. For the case of 2d conformal models it is 
possible to derive exact results. The renormalization group equation 
for the free energy is a reaction-diffusion equation, which is noise- 
free KPZ equation with an additional linear term. For the case of 
Liouville model and strings these models on trees may be naturally 
expressed via the Random Energy Model. This correspondence is used 
to identify the phase structure of strings for analytical continuation of 
DDK expressions. A phase transition is found for spherical strings a 
bit below three dimensions. 



1 Hierarchic tree with branching number close 
to 1. 

One of the most fruitful ideas in physics is the idea of universality. In fact, 
that the only hope due to which rather artificial models of the present the- 
oretical physics can successfully capture the relevant aspects of the nature. 
We believe, that at the critical point the statistical mechanical system omits 
the secondary details. Usually this concerns the Hamiltonian in the d dimen- 
sional Euclidean space. 

It is proposed while keeping the Hamiltonian fixed to simplify the space 
geometry as much as possible retaining two point correlators and three point 
(for isosceles triangles) correlators. If the action of initial theory consists of 
the Laplacian and a potential, our model feels the space dimension through 



the behavior of Green function 



G(x,x') 
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r{x,x'Y ^ 



The total volume is 




(2) 



where L and a are infrared and ultraviolet cutoffs, r{x,x') is the distance. 
The Euclidean geometry has too many constructions. One can rotate a point 
around some center and circumscribe a close circle. Let us now consider some 
metric space with properties: 

A. For every pair of points there is a distance r(x,x'). 

B. There is some measure at every point d/isix) with the total measure 



C.One can construct a quadratic form with corresponding asymptotic (1) for 
the Green function. 

We start out to construct statistical mechanical models on the simplest space, 
that supports points A-C. We hope, that due to the universahty these models 
will acquire some properties of models in d-dimensional space. To realize this 
program we will use some ideas from the theory of Random Energy Model 
(REM) [1-5]. In ref. [5] a relation of 2d quantum Liouville model to REM 
and to the Directed Polymer (DP) on Gayley tree was established. 
Our present analysis shows, that the connection with REM is not a peculiar 
property of Liouville model and works well also for other conformal mod- 
els. Besides, using similar ideas we intend to construct general 2d quantum 
models in the ultramctric space and thereby generalize the above-mentioned 
connection between the quantum field theoretical models and those defined 
on the hierarchical lattices. First, an ultrametric space (surface) with fields, 
located on the surface will be constructed. We define the following three ge- 
ometrical objects: the distance between surface points, the surface measure 
and the volume measure for the ball, delimited by surface. These construc- 
tions are enough to define free field action with correlators that are identical 
to those in the corresponding two-dimensional space. One can add also the 
interaction term to the action. The practical merit of the proposed approach 
is that under certain conditions (for example, Goulomb gas approach to con- 
formal theories) the theories on the ultrametric space can be solved much 
easier, as compared to the Euclidean 2d space. 
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Let us define the ultrametric (UM) space with some measure for the (spher- 
ical) area dns{V,^), volume (ball) measure dny{V,X), total surface and 
total volume — 1. It is commonly known from the mean field theory that 
the surface is of the same order of magnitude as the volume. We can con- 
struct this UM space as a limit of hierarchic lattices. Consider a tree with a 
constant number of branching q in each node and number N for hierarchies. 
The number of end points is , the number of branches - {q^ — l)/{q — 1). 
We consider a set of end points as a surface of sphere, the set of branches 
making a volume (of the ball). Each point on the surface is connected with 
the origin (zero level of the hierarchy) via a single path, which consists of 
links. We determine a measure dfisiV,X) = 1 for each end point, then 
diJ,i{V, I) = {q — 1) for each link. This is mathematically correct in the limit 
N ^ oo. Instead of integer q consider the limit q — > l,A^lng — > Vq. Now 
we have for the total area fj,^ = J (i/i,(Vo,X) = and for the total volume 
Hv = J dfii{V(),l) = e^" — 1. dfiy = J^idfJ'i, where the sum is over the links 
originated from the point X. Later we will manipulate only with dfii. We 
choose g — > 1 to have an equation ni + 1 = /ig. We determine UM distance 
between two points x, y on this surface Vo — v = Vq — jf, where n is the 
number of the hierarchic level, on which x and y had the last common node 
on trajectories to their point from the origin. The maximum UM distance 
between two points v on the surface is Vq ordinary distance as a function 
of V will be defined lately). Now a scalar field (j){x) is defined on our surface. 
For determination of kinetic energy, which should be a quadratic form with 
the Laplacian as the kernel in the conventional space, let us consider the 
expansion 



Here f{v, I) is determined on the links. The integration in (4) is made along 
the trajectory of point X. Since the measures on both the threads of (4) 
coincide (/ dfii(y,l) = 1 + / dVdfii{V,l)), the Jacobian is equal to one. Now 
determine the kinematics part of the action for the field (j){x) 




(3) 




(4) 



Then the partition under the potential U {4>) 



J dfeM-£°dV J dfii{V,l)^f{x,Vf} 
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ex.p{ldf^s{Vo,X)U{cf>{X))} (5) 

We have for the correlator 

< >= V (6) 

where v is the UM distance between the points X, X'. For usual 2d models 
with 

/ #o#exp{- — da;2V0(x)2}exp{ /cia;C/(0(x))} (7) 
J Stt J 

the total surface area is equal to R^, and the correlators read as 

<0(X)0(X')>=ln^ (8) 

In Eq. (7) it is possible to take n-componcnt fields instead of the one- 
component one (f){x). We can determine the distance from the equality 
^ = In r^. Then our correlators coincide (at any rate r ^ 1). It is possible 
to construct a quantum field theory in this case. Our constructions for the 
measure and distance are sufficient. One should bear in mind only, that the 
volume measure inside the sphere V is — 1. 

We are going to discretize (5), then derive the iteration equations for 
imitations of 2d and 3d cases. For the case of spin models on hierarchic 
trees it is well known that it is possible to write simple iteration equations 
(similar to those in Ref.[6]) for any value of q. For the case (5) we formulated 
iteration equations in Section 2 for 2-d imitation and in Section 3 for 3-d case. 
The resulting equation are similar to the KPZ equation [7], but there is an 
additional linear term. Instead of our abstract approach {q = I trees) it is 
possible to consider a branching diffusion process hke [4] and formulate 2-d 
models imitation on its basis. This approach is constructed in Section 4. 

It is possible to derive the majority of results of 2-d conformal theory by 
means of Coulomb gas representation [8]. Here the free field action is modified 
by an imaginary linear term in the action. Correlators of fiuctuating field 
in a critical theories with interaction are equivalent to correlators (from the 
exponential of a such free field) in modified free field picture. In Section 5 
we derived the three point correlators in our approach and the results agree 
with those in the 2-d case [8]. 
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For the case of exponential function for U (x) potential (2d Liouville model 
in the Euclidean space) there is a strict result [9]-[10] that the thermodynam- 
ics of our model on hierarchic tree is independent of q and is similar to REM. 
In Section 6 a qualitative derivation for a REM with complex replica numbers 
(using the results of [11] for real number of replicas)is given and the solu- 
tion of REM at complex temperatures [12]- [13] is obtained then those results 
are used to identify phase structure of strings, using analytical continuation 
of DDK formulas [14]- [16]. In Section 6 the main results of the work are 
discussed. 



2 Iteration equations for the 2-d case. 

The advantage of representation (5) is that we are in a position to calculate 
the partition function through iterations. This is well known for models on 
hierarchical lattices [6]. Let us for some large number K, divide V into K 
parts V/K and determine a hierarchical tree with K levels and branching 
number 

? = exp{^} (9) 

Similar to [5] it is easy to define the partition Z via iterations. In case of 
some large number K we derive 

/ K 1 
= J_^exp{--Ky^ + U{x + y)}dy 

li+iix) = y — J_^exp{--Ky'}[I,{x + y)rdy 

Z = lim [Ik{0)Y (10) 

K—^0O 

As for the determination of partition function, we need only the equation 
(10). Our choice g ^ 1 is a reasonable simphfication. It is possible to 
construct pcrturbativc field theory, calculate diagrams. To solve analytically 
equation (10) it is convenient to consider the other (opposite) case, q oo. 
Only in this case the analytical solution proves possible. The point is, that 
for the bulk structure of theory the value of q is irrelevant. For 
example, in asymptotic expansion of free energy 

F(iV) = FoiV + FilniV + F2... (11) 
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only the last term depends on the choice of q (in case of Directed Polymer 
N ^V). 

We found how to construct the simplified version of any 2d theory on a 
hierarchical lattice. 

Let us consider carefully equation (10) in the limit of large K and in- 
troduce variable w{v,x) = Ikv (x). We consider the limit ^ <^ 1. For the 

differential dv we have an expression ^. Let us also assume 

q-l = ^ = dv (12) 

Using the expression f» x{l + logx{q — 1)) it is easy to obtain 

dw 1 
—— — winw + - Aw 
dv 2 

w{0,x)^exp{-U{x)) (13) 
After the replacement w = exp{u{t,x) we arrive at 

u{Q,x) = U{x) (14) 

where U{x) is the potential in Eq.(7). Having an expression for u{v,x), we 
obtain for the free energy 

In Z = u{V,0) (15) 

We have a noise- free KPZ equation (14) with additional linear term for the 
free energy. 

There are two interesting solution of Eq. (14) at large values v. If the 
couplings in the polynomial potential are 0(l),it is reasonable at large values 
of V to consider the solution: 

u{v,x) = const exp{v) (16) 

If one considers the couplings ~ gyp^^y^ in the potential U (x) , then the 
solution 

u{v,x) = const exp{v) +Us{x),Us{x) ~ 1. (17) 
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^Aus + ^{Vusf + Us = (18) 

corresponds to the perturbative regime. This equation gives the effective 
potential at the stable point of renormalization group. One can rewrite of 
Eq. (18) in another form for z = 

dz 2 

— + Z + -U, 19 
dUg z 

In analogy to Eqs. (10), (13) it is possible also to derive the correlators. 
To calculate the correlator < exp{ia(l){x) — ia((){y)) > , where the hierarchic 
distance between points x, y is vq, one should distinguish during the iteration 
between the links located on the paths that connect the origin with the points 
x,y. Thus we also consider the equation 

^ f\nw + -A/ 

dv 2 

/(O, X, a) = exp{U{x) + lax) (20) 

Then for the generating function /o(f,x) of correlator one must solve again 
Eq. (20) with the boundary conditions at the point vq 

fo{vo, x) = f{vo, X, a)f{vo, x, -a)/w{vo, x) (21) 

We obtain an expression for the correlator : 

< expiiacpix) act>{y)) >= (22) 

In this way we can calculate two point correlators, as well as other multipoint 
correlators. 



3 High dimensions 

The same approach may be used for the case of d > 2. In d-d space one has 
for the volume ~ a''. If we identify it with our q^, then a — qd^. To have 
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exact expression for the correlator, // are defined on the branches, /o, /i at 
the origin. Here the free field action is defined as 

0(^) = /o + /i + E/KO (23) 

vl 

The summation in (23) is along the trajectory of point X. Now determine 
the kinematical part of the action for 0(a;) field 

^ = ^[/i' + E exp(-a^)/(t;, if /a] (24) 

^ vl 

If one takes a — ^ for the combined field, then 

< 0(x)0(x') >= exp(Q;^;) ~ [— ^]-(<^-2) (25) 

r(a;, x') 

where L is the infrared cutoff. Now (10) transforms into: 

^ = ii'oo V 2^ i-J''^^~2y' + U{y)}dy[lK{y)n2Q) 

To calculate Ik{x) = w{V,x) = exp{u{V,x)) we have to solve the equation 
hke (14) 

^ = \d{d-2)exp[d{d-2)v]Au+\{Vuf + u 
dv 2 2 

u{0, x) = U{x) (27) 

It is important to investigate the version of Generalized Random Energy 
Model (GREM) corresponding to (27) , when our model is defined on hier- 
archic lattice with large branching number q. Here the physics at d = 2 and 
d > 2 is quite different and one cannot use the methods of [9] for the latter 
case. 
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4 Branching diffusion 



Instead of the structure oiq — 1 trees it is possible to consider an ensemble of 
hierarchic trees having again a small branching number after the averaging. 
An analogous process has been introduced in [4], and here we give only a little 
modification to choose proper boundary conditions for a desired potential. 
There is an origin and a branch from it. Branch appears from the original 
branchs during the period of time dt with the length dt. Alternatively, with 
the probability 1-dt the old branch is elongated by the length dt and a random 
variable // is introduced with the variance 

< /f >= dt/2 (28) 

After some period of time t the number of endpoints is exp(t). Every branch 
has one or more random variables. The fields at the endpoints of tree are 
again defined as a sum of random variables along the trajectory. The variable 
w{t,x) is defined as 



w 



{t,x)^<exp{J2U{y + x))> (29) 



Here the summation y is along all endpoints of a tree after period of time 
t. Let us define w(t + dt,x). One should keep in view at the determination 
of w{t + dt,x) that during dt time either there appears a new branch (the 
contribution of this process is equal to dtw{t, x)"^), or , if stays the old branch 
has an increment of the contribution is {1 — dt)w(t,x + fi). Combining 
the contributions we easily derive the KPP equation [9] for a : 

^la^ = w^^t, x) - w{t, x) + -Aw (30) 
dt 2 

The derivation of this equation has been done in [9] where the case w(0, x) — 
exp(— e^) has been considered. It is similar to (13), with the only difference 
that the nonlinear term is replaced by — w. We believe that critical prop- 
erties of both the approaches (13), (30) are the same. The point is that for 
q = 1 tree approach there is no any averaging in an ensemble and thus it is 
easier to deal with the perturbative expansion. Of course our equation (13) 
is also a reaction diffusion-equation. 

Another interesting problem is an imitation of d-d reaction diffusion equa- 
tions on q = 1 space. 
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5 Coulomb gas representation for 2d confor- 
mal fields. 



We can apply these ideas to conformal theories using the Coulomb-gas for- 
malism with the background charge ao [8]- If we correctly defined the zero 
mode of Laplacian and the correlator has correct dimension, we have good 
chances to imitate the 2d situation. We have an action 

[ d^wJ§(l)A(j) + i2V2aoR(p (31) 

Here the field (^{w) is defined on the sphere, R is curvature, A is a Laplacian. 
One defines the screening charges from the condition a'j_^ — 2aQa±i = 1. To 
calculate the correlator Hit exp{i\/2Q!fe0} with screening charges Q'^Q'^ one 
has to consider 

Z ^ J Dg(f)e^ -f 'i'^'>^V^'P^<t>+i'^V2aoR4> 

llcxp{iV2ak(p}i J d^w exp{iV2a+(py\ J d'^wexp{iV2a_<f))"' (32) 

The zero mode integration gives for nonzero correlator the constraint 

^ Q!j -|- ma^ + na^ — 2ao (33) 

i 

in 2d case [8]. 

To have finite set of ctj from (35) we should put a constraint 

p'a+ + pa^ = (34) 

which is the definition of minimal models. The deficiency of our approach is 
that we cannot find connection between ao and conformal charge c. 
While calculating (26) in UM space, wc omit the ao term. We again , as in 
2d case, consider normal ordered operator product for exp[i 

< z'^{V2a+)z''{V2a_) exp[iV2aEi] exp[iV2aE2] > (35) 

Here the average is over the normal distribution on our hierarchic tree. As 
m and n are integers, it is possible to perform the integration via directly. 
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How one can derive the expression for the pair correlator? All the n + m 
charges with the total charge 2ao — a are located near the point 1, or the 
point 2, so we have for (35) 

< z'^{V2a+)z''{V2a-)exp[iV2aEi]exp[iV2aE2] >~ 

exp[a{2ao — a)v] — exp[a{2ao — a) In — ] (36) 



For 3 point correlators we now consider expressions like 

< z"'{V2a+)z''{V2a^) e:K.p[iV2aiEi] exp[iv/2 0:2-^2] ex.p[iV2a3Es\ > 

vn = V23 > vl2 = V (37) 

We assume that the screening charges are near the points 1 and 2. Then the 
distance from the point 3 to any charge is equal to R and we immediately 
obtain for the dependence of (37) on R: 

^-2a3(2ao-a3) ^3g^ 

Let us consider the dependence of (32) on r = Lexp[|]. 
As a result of direct interaction of charges ai, a2 we obtain: 

If all our n + m charges are located at the distance r from both the points 
1,2, then their interaction energy with charges 1, 2 is proportional to 

2{ai + a2){na+ + ma^) (40) 

For the energy of self-interaction of the screening charges we have: 

— {na^ + ma-Y — {n-\-m) — 2aQ{na^ + maJ)] (41) 

The integration over the coordinates of n + m charges gives exppv], so even- 
tually we have: 

^— a3(2ao— a3)j,2aia2 + (na++ma-)^— 2ao(na++ma_) (42) 
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It is easy to check that this expression is equivalent to standard expression 
from the conformal field theory 

^-2a3{2ao-a3)^al-2aiao+al-2a2ao-{ — 0:3+20300 (43) 

Combining two expressions, we derive eventually for the correlator: 

203(200— as) j,2aia2+(na++ma-)^—2ao(na++ma_) (44) 

Let us put in [43] ru — ri3 = r23 in (43) and consider REM instead of the 
Directed Polymer. In case of REM we have that < (f){X)(f){X') >= 5{x - x'). 
It is possible to investigate the phase structure of (originally DP) correlator 
in this way. 



6 REM at complex temperatures with com- 
plex numbers of replicas and strings. 

6.1 REM version of strings 

Using similar ideas we will connect a string partition (after integration by 
zero mode) with finite replica REM and investigate the phase structure. 
Recall some results from string theory for a string in = c space with 
spherical surfaces [14-16]. It is known for the partition that after integration 
by zero mode 



/ 



where 



c 



- = ^[c-25-V(25-c)(l-c)] 



(46) 



Here 4>{w) is a field on two dimensional sphere with coordinates w, curvature 
R. If we continue those formulas for c > 1, the coefficients become complex. 
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We see, that Z ^< z'^ >, where the averaging is over the normal distribution 
of field (p{w). 

If the conjecture about the equivalence of 2-d model (49) and corresponding 
q = 1 model is correct (it could be checked numerically), then one can solve 
explicitly DP problem in Q ^ oo limit [9]- [10]. In this limit it is easy to 
prove [10], that the Directed Polymer has the same thermodynamic limit 
(Fq), as a simple REM. 

If in (45) our variables (f){w) are distributed according to the normal law 
J d^w^g<t>A4>+QR4> ^j^-^ non-diagonal quadratic form, in the case of REM 
all (piw) are independent variables with a normal distribution. If we replace 
the model (45), (46) with models on g = 1 trees, then according to [9]- [10] 
these models arc equivalent to REM. Thus it is worthwhile to consider REM 
instead of the set (45)-(46). 

To construct an equivalent scheme of REM let us introduce infrared and 
ultraviolet cutoffs L and a. Then the physical number of degrees is 

7-2 

M = — (47) 

Now define the distribution of 4>{w) over all points w, using the free field 
action from (45): 

p(0o) =< 5(00 - cl){w) >^(^)~ exp(-^^) (48) 
where the averaging is over the distribution 

p(0(w)) - 68^ / d^w^9<l>A^+QR^ (49) 

and 

G'(0) = 21n- (50) 
a 

We can replace our system with a collection of M independent variables 
Ei ~ 0(w) with the distribution (43) instead of (44). Our goal is to calculate 
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It is possible to solve the system rigorously. Here we are giving a qualitative 
derivation, (checked by our exact calculations). 
Note that 



6.2 Solution of REM at complex replica numbers 

Let us consider Eq. (51) for positive integer values of /i, where the averaging 
is made over the distribution (52) for each Ei. 

There are only two competing terms in the sum and two corresponding 
phases. The first one is paramagnetic (PM) phase , originated from the 
cross terms in expansion 

lnZ = ^lnM + A'^ = ivM±i3ii (53) 

The second one, the correlated paramagnetic (CPM)[11] is originated from 
the diagonal terms in (53) like e"'^'^^' 

M 

Z=< (^e-'^'^^O > 
lnZ^lnM+^^^k±dV) (54) 

Let us consider the continuation of (54) at < 1. At critical /3c its entropy 
InZ — P^^^ disappears. Let as assume for InZ in this region an expression 
proportional to (3 (it is natural for a system with zero entropy) and /i. The 
continuity of In Z gives for the spin-glass (SG) phase 

In Z = NiiPcP (55) 

If we pass to complex temperatures [14]-[15], then (51) transforms to (it is 
easy to check this directly for integer /x) 
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In Eq. (55) one has to replace /3 by /3i,then 

In Z ^ Nii(3^(3i (57) 

For complex temperatures there is the fourth, Lee- Young-Fisher (LYF) phase 
[12]. Its derivation is not direct. The point is that for noninteger values of fi 

Z^<\z\''> (58) 

After this trick it is easy to derive the CPM expression. The principal terms 
are now given by terms 6"^^^^* 

l„^ = iM±M)^^ (59) 

4 

Now continue our four expressions to complex values of /i. 
For PM phase 

,„^^^M±iiziSK (60) 



For SG phase 



For LYF phase 



For CPM 



\nZ ^ NfjLiP^Pi (61) 



,^,^m±m^ (63) 

The imaginary parts of InZ in (61)- (63) were ignored. 

To find the borders between four phases first has to obtain the correct phase 
at ^ limit, then compare its expression for | lnZ| with the corresponding 
one for CPM phase. It is known that LYF phase exists at 

/5 < I (64) 
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and PM one at /? < (3c- For complex temperatures one has a condition for 
SG phase 

(3i>(3c + P2 (65) 

The last point. The rigorous derivation gives that LYF for noninteger //i 
exists only at 

Hi > -2 (66) 
6.3 REM results for the phase structure of string. 

One can apply this (though qualitative but rather strict) result to strings. 
Identifying ^ P,^ ^ 1-^ and using (2) we derive 

/3i = P2 = Hi = ^[25-c],H2 = V(25-c)(c-l)^ (67) 



For other string topologies one has 

H^{i-g)H (68) 

Note that y = LYF phase exists only in the torus case at 25 > c > 19. 
At 19 > c > 1 the system with torus topology and higher is in SG phase. 
For sphere topology 19 > c > 1, it is in CPM at 

l + 4y(2y-l)2>4yi (69) 

otherwise in SG phase. We see phase transition at d ~ 2.98. 
For 25 > c > 19 spherical topology case CPM is at 

l + 4y(2y-l)2> i + V, (70) 

otherwise at LYF phase. What can we say about string physics base on the 
REM picture? 

The most interesting case is the spherical one. When the value of as in- 
creased to pass over the c = 1 barrier, nothing happens in REM picture, the 
system is still in CPM phase, as for in the c < 1 case. The PM or CPM 
phases are ordinary physical phases, so one could try to succeed here with 
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the same level of reliability, as for c < 1 strings outside of the minimal se- 
ries. To reveal interesting (unitary) theories explicitly one should solve the 
directed polymer at finite replica number including finite size corrections and 
correlators. 

For the high topologies we have omitted modular space dependence of par- 
tition. It will be interesting to check our conclusion about different physical 
phases for different topologies numerically. But at least for spherical case the 
REM analysis seems quite rehable. 

7 Conclusions 

We generalized the results of [5] , obtained for the 2d model with an exponen- 
tial potential, for the case of other critical models, that may be formulated by 
means of Laplasian (as a single differential operator) and a desired potential 
of any dimensions d>2. The models can be formulated on hierarchic trees 
with constant branching number q. Main observation of this work is that 
for the special case of g — 1 it is possible to construct a field theory that is 
similar to renormalized field theories in continuous spaces. 

We hope, that the bulk structure, the two and three point correlators 
(for isosceles triangles) are the same, as those in d-d critical models. This 
hypothesis is correct for the case of Liouville 2d model , as well as for the 
(free) case of Coulomb gas representation [8] . 

We have also analyzed the phase structure of strings by consideration 
of an analytical continuation of DDK formulas by analogy with REM and 
solution of complex replica REM at complex temperatures. 

It is possible to check our hypothesis about the equivalence of our models 
on q — 1 trees with some segment of d-d field theory by means of direct 
numerical calculation of Eq. (14), (27), for example for, the field version of 
3d Ising model with proper choice of potential U. 

I am grateful to ISTC grant A-102 for partial financial support, C. Lang and 
W. Janke for invitation to Graz and Leipzig and discussions, P. Grassberger 
for useful remark. 
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